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ABSTRACT 

The problem of relativistic stellar pulsations is studied in a somewhat ad hoc approx- 
imation that ignores all fluid motions. This Inverse Cowling Approximation (1CA) is 
motivated by two observations: f ) For highly damped (w-mode) oscillations the fluid 
plays very little role. 2) If the fluid motion is neglected the problem for polar oscillation 
modes becomes similar to that for axial modes. Using the ICA we find a polar mode 
spectrum that has all features of the w-mode spectrum of the full problem. Moreover, 
in the limit of superdense stars, we find the ICA spectrum to be qualitatively simi- 
lar to that of the axial modes. These results clearly show the importance of general 
relativity for the pulsation modes of compact stars, and that there are modes whose 
existence do not depend on motions of the fluid at all; pure "spacetime" modes. 
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1 INTRODUCTION 

The pulsations of relativistic stars have attained increasing 
interest in the scientific community in the last four years. Re- 
cent results have improved our understanding of both this 
ject, and the theory of general relativity itself Never- 



liable results and also reveals the limitations of the specific 
description of the pul sation problem that has been used i n 
all calculations so far (Andersson, Kokkotas & Schutz 1995). 



The oscillations of relativistic stars are often described 



subjc 

t heless , I one can easily identify issues that are far from well 
understood at the present time. 

When Thorne and his colleagues established th e theory 
for relativistic stellar pulsations three decades ago (Thorne 



by a system of four coupled ODEs ( Detweiler & Lindblom 



|fc Campolattaro 1967[ [Thorne 1969| ), the main results were 



not that different from the well-known ones of Newtonian 
theory. Compact objects where found to oscillate at almost 
exactly the frequencies that Newtonian theory predicted. 
Relativistic effects only gave rise to a very slow damping of 
the pulsation. The emission of gravitational radiation im- 
plied that the oscillation frequencies were complex with a 
relatively tiny imaginary part. 

The situation recently changed considerably when a 
new family of oscillation modes for comp act stars was found 



1985). Two equations correspond mainly to the fluid pulsa- 
tions and the other two to perturbations of the spacetime. 
The two sets of equations couple strongly in the case of 
highly relativistic stars, but the coupling also depends on 
the frequency. For slowly damped modes, the equations are 
in practice decoupled in the high and low-frequency lim- 
its. These limits correspond to the fluid p and (/-modes, re- 
spectively. The two spacetime ODEs then hardly affect the 
structure of the spectrum at all, and a relativistic g eneralisa- 
tion of the Cowling approximation (fowling 1941), in which 



perturbations of the spacetime itself are igno red, yields the 
correct spectrum [see for ex ample Robe ( 1968 ), McDermott, 
Van Hor n and Schorl ( |l983| ), Finn ( |l98cj ) and Lindblom and 
Splinter (|l990|)]. 



(Kojima 1988; Kokkotas & Schutz 1992). These modes have 



no relation to the known p and (/-modes in Newtonian the- 
ory. Rather, this new family of modes, the exist ence of which 
was suggested by a simple model problem (Kokkotas & 



([l99§) 



Inversely, as has been shown b y Kok kotas and Schutz 
Leins, Nollert and Soffel ( 199J ) and Andersson, 
the fluid hardly pulsates at all 



Kokkotas and Schutz (199 



Schutz 1986), comes from coupling the stellar fluid to the 



spacetime. A characteristic property of the new family is 
that the fluid motion is very small, and the oscillations are 
rapidly damped. These modes have been named w-modes 
(gravitational wave modes), and their existence (together 
with an additional part of th e spe ctrum) has been verified 
by Leins, Nollert and Soffel (1993). Moreover, we have re- 
cently developed an accurate numerical code that gives re- 



in the case of ui-modes. A simple way to find the spectrum 
of such modes might therefore be to omit the equations per- 
taining to the fluid perturbations altogether; what we will 
from now on refer to as the inverse Cowling Approximation 
(ICA). In the present work, which is intended as a comple- 
ment to studies of the full problem, we examine the stellar 
pulsation problem in this approximation. The major motiva- 
tion is to clearly show the degree of involvement of general 
relativity. The ICA should indicate which features of the 
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w-mode spectrum that are related to the fluid motion and 
which are due to the spacetime perturbations. This is an im- 
portant step towards a better understanding of the nature 
of the w-modes and the properties of rapidly damped stellar 
oscillations. 



equations imply the relation ( Detweiler & Lindblom 198 



Chandrasekhar and Ferrari (1991) recently showed that 
axial quasinormal modes could exist for very compact stars. 
The possibility of axial modes had previously been discarded 
since axial perturbations do not couple to oscillations of the 
stellar fluid. However, Chandrasekhar and Ferrari realized 
that "trapped" modes can occur when the star is so compact 
that the surface lies inside the peak of the famil iar Regge- 
Wheeler potential barrier ( Chandrasekhar 1983 ). That is, 
the axial modes can be understood in terms of a po tential 
well [see figure 1 in ( Chandrasekhar fc Ferrari 1991 )]. Ini- 
tially, only a few such axial modes were found and t hey were 
all slowly damped. Recent work by Kokkotas ( 1994 ) revealed 
that there are a large number of highly damped axial modes 
as well. This new set of axial modes is in many ways similar 
to the polar w-modes. 

A good reason for studying the axial modes is that the 
spacetime perturbations do not couple to those of the fluid, 
and thus the system of equations becomes a very simple one. 
In fact, the problem is similar to that for a perturbed black- 
hole, although the boundary conditions (now at r — » 0) are, 
of course, different. If we freeze the fluid motion in the case of 
polar perturbations we arrive at a similar problem. Using the 
Inverse Cowling Approximation we can consequently study 
the polar pulsations in a way similar to that used for the 
axial ones. This should provide a better understanding of 



the relation between the two sets of pulsation modes for 



relativistic stars. 

Recently, two of us have been involved in a study of 
superdense stars cl ose to the limit of compactness posed b y 
general relativity (Kojima, Andersson & Kokkotas 1995). 



That study suggested that the origin of the axial modes and 
the polar tii-modes is the same. Both sets are "spacetime" 
modes that do not depend on the stellar fluid at all for their 
existence. If that conclusion is correct one would expect all 
features of the ui-modes to be present in the ICA. 



2 THE INVERSE COWLING 
APPROXIMATION 

In Regge- Wheeler gauge the perturbed metric can be written 

ds 2 = — e" (l + r Hoe^Yim) dt 2 — 

2iuir +1 H\e Mt Yimdtdr + 

e x (l - r H e wt Yi m ) dr 2 + 

r 2 (l - r e Ke'^Y em ) (d6 2 + sin 2 6d<t> 2 ) , (1) 



with Ho, Hi and K functions of r only. Ye m are the stan- 
dard spherical harmonics and M(r) acts as an effective mass 
inside radius r. It is assumed that all perturbations depend 
on time as exp(iajt). We also have 

.--i-^. P) 

The metric function u, the pressure p and the mass M fol- 
low from the Tolman-Oppenheimer-Volkov equations that 
determine a stellar equilibrium model. Moreover, Einstein's 
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H 
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Hence, only two metric functions remain undetermined. 

The full pulsation problem consists of four coupled first 
order ODEs [equations (8)-(ll) of Detweiler and Lindblom 
( 1985|) 1 These describe the two remaining spacetime vari- 
ables Hi and K as well as two fluid ones: V and W. In this 
formulation, the ICA corresponds to V = W — and the 
perturbation equations for a specific multipole £ take the 
following simple form inside the star; 



rH{ = e A [Ho + K] 
rK' = H +£(£ + i)Hi 



1 H e + 4nr e (p — p) 

r 



£+1 



[M + 4vrr 3 p) 



Hi ,(4) 
K,(5) 



where a prime denotes a derivative with respect to r and (3) 
should be used to replace Ho- 

In the exterior the perturbati on equations s i mplify 



to the well-kno wn Zerilli equation (Fackerell 1971, Chan- 
drasekhar 198J). It is, in fact, possible to transform equa- 



tions (4) and (5) into a single second-order equation but we 
have found the form of this equation to be rather compli- 
cated and not numerically convenient. Instead, we approach 
the above system of equations numerically in the present 
study. The procedure used is a very simple one: The result 
of one integration [initiated with the regular solution to (4) 
and (5)] from the centre of the star to the surface is matched 
to an exterior solution calculated according to the m ethod 
described by Andersson, Kokkotas and Schutz (1995). The 



technique used for finding eigenfrequencies here is identical 
to that discussed in our previous paper. 

Before we discuss the numerical results obtained in this 
way a few words of caution are in order. We have not de- 
fined the ICA in a gauge-invariant way here. Such a defi- 
nition does, in fact, not seem possible. This does not mean 
that the idea presented here is without merits, however. The 
usefulness of an ad hoc approximation such as the ICA is 
illustrated by the actual results obtained. Even though such 
results should not be used as indication of new physical phe- 
nomena we will show in the following sections that they pro- 
vide interesting information. We view the present study as a 
mathematical experiment that provides information on the 
relative importance of the spacetime variables and the fluid 
ones. 



3 ICA MODES FOR POLYTROPES 

We have calculated the quadrupole quasinormal-mode spec- 
trum in the IC A for the four stellar models of Kokkotas and 
Schutz (11993). The results are interesting. First of all, we 
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Figure 1. The ui-mode spectrum compared to the ICA spectrum 
for a simple polytropic stellar model. Here, the equation of state 
is p = 100 km 2 p 2 , and the star has characteristics: p c = 3 X 10 15 
g/cm 3 , R = 8.861 km, M = 1.266M Q , i.e., 2M/R = 0.422. This 
corresponds to model 2 of Kokkotas and Schutz (1992). 



could not find the p and (/-mode spectra. This was certainly 
expected since these modes are clearly associated with fluid 



pertur bations which are absent in our approximation ( Mc- 
Dermol|t, Van Horn fc Scholll983|) . 



In Figure 1 we compare the ui-mode spectrum to the 
ICA mode spectrum for model 2 of Kokkotas and Schutz 

(1992) . The similarities between the two spectra are evi- 
dent. The highly damped part of the ICA spectrum, i.e., 
the w-modes, is in excellent qualitative agreement with the 
results of an analysis of the full problem. In fact, the real 
parts of the frequencies agree surprisingly well. The spacing 
between consecutive modes is roughly 10 % smaller in the 
ICA than in the full problem. This means that the abso- 
lute difference between the real parts for "corresponding" 
modes in the ICA and the full problem increases drasti- 
cally with the frequency. Nevertheless, the relative difference 
[(Re u>n+i — Re oj n )/Re lu„+i where n is an integer labelling 
the modes] is typically smaller than 10 % . As for the damp- 
ing rates, the ICA imaginary parts are some 20-30 % smaller 
than those for the w;-modes of the full problem. As the oscil- 
lation frequency increases this difference decreases and the 
two spectra approach each other. All characteristic features 
of the spectrum, such as the existence of a few overdamped 
modes with small real part [what Leins, Nollert and Soffel 

(1993) refered to as wrr-modes] are found where expected. 



Even though a discussion of physical effects here must 
carry a disclaimer because of the gauge-dependency of the 
ICA it is interesting to speculate on reasons for why the 
damping of the ICA modes is slower than that of the w- 
modes of the full problem. A hand-waving argument is the 
following: In the full problem the fluid is also "radiating" 
gravitational waves with the signature of a ui-mode (the 
perturbations are coupled). The fluid thus acts as a kind 
of "gravitational-wave pump" and it seems plausible that 
the dissipation rate increases. The observed result would fol- 
low if the fluid is "a more efficient radiator of gravitational 
waves" than the spacetime itself. But does not the evidence 
from oscillations mainly associated with the fluid (the p- 
modes) indicate the opposite? The damping of these modes 
due to gravitational radiation is slow, and the fluid would 
seem to be a poor radiator. This is certainly true, but it is 
important to remember that the origin of the fluid modes 
and the "spacetime" modes we discuss in this paper is quite 
different. For the fluid modes the spacetime perturbations 
are negligible compared to the fluid ones, and the modes are 
slowly damped because the coupling between matter and 
gravitational waves is weak. The situation is different for 
the w-modes, for which the spacetime perturbations play 
the dominant role. The results in Figure 1 suggest that the 
spacetime curvature traps gravitational waves more effec- 
tively than does the stellar fluid. This seems plausible, but 
further study of this problem is needed if we are to under- 
stand the actual physics involved. 



4 ICA MODES FOR SUPERDENSE STARS 



Let us now compare the polar modes obtained in the ICA 
to the axial modes. To do this we consider the uniform den- 
s ity m odel that w as stu died by Chandrasekhar and Ferrari 
( 199l j ), K okkotas (1994) and Kojima, Andersson and Kokko- 
tas (1995). For this model the surface of the star can be 
inside the peak of the curvature potential. In this way we 
have a problem that can loosely be characterized as a "po- 
tential with a well inside a barrier", and one would expect 
modes with small imaginary part (which come mainly from 
the potential well, i.e., act as quasi-bound states in quantum 
l angua ge) to exist. Furthermore, as was shown by Kokkotas 
(1994), there are modes with high damping. 

Before we proceed to discuss the present results it is 
necessary to discuss the discrepancy between pr eviou s re- 
sults for axial modes. In the study of Kokkotas ( 1994 ) the 
damping rate of the modes was found to be r oughly half 
that found by Chandrasekhar and Ferrari (1991). When in- 
vestigating this issue we found it to be due to a misp rint 
in equation (19) of Chandrasekhar and Ferrari (1991). If 
the erroneous equation is used in the numerical calcula- 
tions the results of Kokkotas follow. Once the equation is 
corrected the numerical results are in good agreement with 
those of Chandrasekhar and Ferrari. It is, however, impor- 
tant to stress that the results of Kokkotas are qualitatively 
correct. An infinite spectrum of highly damped modes does, 
indeed, exist. 



Here we have used the method of Andersson, Kokkotas 



and Schutz (1995) to approach the exterior problem and ver- 
ify the qualit ative b ehaviour of the recent axial-mode results 
of Kokkotas ( 1994 ). It is worth noticing at this point that, 
altho ugh t he numerical results of Chandrasekhar and Fer- 
rari (1991) are correct their method of outward numerical 
integration is very sensitive to the choice of endpoint rep- 
resenting "infinity" . Integration from the surface of the star 
towards infinity is, in fact, not very reliable in this kind of 
problem, especially not when one deals with numbers of the 
order of 10 -8 as one must to identify the long-lived modes. 
Integration towards the surface of the star , as in the method 
used here, is considerably more stable and accurate. 

For these ultracompact models we find that the ICA 
spectrum is quite si milar to the axial-mode spectrum ob- 
tained by Kokkotas (1994), see Figure 2. The ICA modes 
are generally slower damped than the axial modes through- 
out the spectrum. This is reminiscent of the result for po- 
lar modes of polytropes discussed in the previous section. 
Our calculations have also unveiled highly damped modes 
with relatively small real parts in both spectra, see Figure 
2. These modes are, in many ways, similar to the "new " po- 
lar modes identified by Leins, No llert a nd Soffel (1993) and 
Andersson, Kokkotas and Schutz ( 1995| ). That similar modes 
exist also for axial perturbations was not known previously, 
but is evident from our Figure 2. 

The slowest damped modes for very compact stars 
should be view ed as modes trapped inside the curvatu re po- 
tential barrier (Kojima, Andersson & Kokkotas 1995). The 



slight difference in damping between the ICA modes and the 
axial modes in the first part of the spectrum could be due 
to the difference between the Regge- Wheeler and the Zerilli 
potential. Although similar, these two potentials are not the 
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Figure 2. The axial and ICA spectra for a very compact uniform 
density star. This specific example is a star for which R/M = 2.28 
i.e., 2M/R = 0.88. The modes that correspond to the "new" polar 
modes identified by Leins, Nollert and Soffel (1993) lie above the 
general "string" of modes. That such axial modes exist was not 
known previously. 



same ( Chandrasckhar 1983) and a difference in the stellar 
spectra should be expected at some level. 

The rapidly damped modes cannot easily be viewed as 
trapped modes in this sense. Rather, they are analogous to 
the w-modes for less compact stellar models discussed in the 
previous section. In fact, it is worth stressing the consider- 
able qualitative similarity between the rapidly damped axial 
modes and the polar to-modes. (compare Figs. 1 and 2). It 
should also be rememb ered that the results of Kojima, An- 
dersson and Kokkotas ( 1995 ) show that the axial and the po- 



lar spectra approach each other as the star gets increasingly 
compact. By extending that study to less compact stellar 
models one may hope to shed further light on the relation- 
ship between the axial modes and the polar ui-modes. Such 



work is presently in progress ( Andersson, Kojima fc Kokko 
tas 199™. At present it seems clear that, alth ough axial 



pertur bations do no t couple to the stellar fluid ( Thorne & 
Campo attaro 1967 ), rapidly damped axial modes should ex- 
ist also for less compact stellar models. That is, models for 
which the surface of the star lies well outside the peak of the 
Regge- Wheeler potential barrier should support a branch of 
strongly damped axial modes. 



5 CONCLUDING DISCUSSION 

In this short paper we have presented results for a new ap- 
proximation relevant to stellar pulsation problems in gen- 
eral relativity. The Inverse Cowling Approximation, which 
neglects perturbations of the stellar fluid, provides a useful 
tool for probing the role of the spacetime degrees of free- 
dom in this problem. We have compared 1) the spectrum of 
highy damped modes for polar perturbations (the «;-modes) 
to the corresponding ICA modes for polytropic neutron star 
models and 2) the polar ICA modes to the modes for axial 
perturbations of extremely compact uniform density stars. 
The results are unequivocal: All essential features of the w- 
mode spectra are present in the ICA. This is strong support 
for the idea that the iu-modes are "spacetime" modes that 
not rely on the motions of the fluid for their existence. 

Motivated by the present results, we would argue that 
this kind of approximation can be used to further improve 
our present understanding of the stellar pulsation problem. 
In fact, it works much in the same way as the Cowling ap- 
proximation does for pulsations mainly associated with the 
fluid degrees of freedom. In that case, calculations are facili- 



tated by neglecting the perturbations of spacetime (McDer- 



mott, Van Horn & Scholll983). When relativistic stars are 



considered, that approximation on its own does not make 
much sense. It is clear, eg. from the existence of w-modes, 
that the gravitational field must be considered as dynamic 
if all features of the problem are to be accounted for. So 
only if the Cowling approximation is complemented with 
something like the present approximation can one infer the 
relevant physics. 
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